Enhanced photon-assisted spin transport in a quantum dot attached to ferromagnetic 

leads 
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We investigate real-time dynamics of spin-polarized current in a quantum dot coupled to fer- 
romagnetic leads in both parallel and antiparallel alignments. While an external bias voltage is 
taken constant in time, a gate terminal, capacitively coupled to the quantum dot, introduces a 
periodic modulation of the dot level. Using non equilibrium Green's function technique we find that 
spin polarized electrons can tunnel through the system via additional photon-assisted transmission 
channels. Owing to a Zeeman splitting of the dot level, it is possible to select a particular spin com- 
ponent to be photon-transferred from the left to the right terminal, with spin dependent current 
peaks arising at different gate frequencies. The ferromagnetic electrodes enhance or suppress the 
spin transport depending upon the leads magnetization alignment. The tunnel magnetoresistance 
also attains negative values due to a photon-assisted inversion of the spin-valve effect. 



I. INTRODUCTION 

Time-dependent transport in quantum dot system 
(QDs) has received significant attention due to a vari- 
ety of new quantum physical phenomena emerging in 
transient time scaleii A few examples encompass charge 
pump^i^ and photon- assisted tunneling transport.—"— For 
instance, a double dot junction sandwiched by leads can 
be used to pump electrons uphill from a lead with lower 
chemical potential to a lead with higher chemical poten- 
tial, in contradiction to the usual dc-regimci^ This was 
achieved by applying a sinusoidal gate voltage among 
the dots. Photon-assisted tunneling can occur when an 
oscillating gate potential or laser field is applied in a 
QD or a metallic central island coupled to source and 
drain terminals4r"— Time-dependent regime also leads 
to zero-bias charge or spin pumping when a minimum 
set of two parameters of the system (e.g., gate poten- 
tial and tunneling rate in a QD system) are time modu- 
lated independently. This is the case, for instance, of the 
non-adiabatic charge and spin pumping through interact- 
ing quantum dotsii and quantum pumping in graphene- 
based structures i^^rJ^ 

Transient charge and spin dynamics in an interacting 
QD driven by step pulse or sinusoidal gate voltages re- 
vealed distinct charge and spin relaxation timesj^ An 
exquisite behavior that has been predicted theoretically 
is the self-sustained current oscillations in a quantum dot 
system driven out-of-equilibrium by a fast switching on 
of the bias voltage, contrasting to the expected steady 
state behavior. This phenomena has been attributed to 
dynamical Coulomb blockade.— 

It is in the fascinating area of spintronics^^ that time- 
dependent quantum transport reveals its prolific poten- 
tiality in producing spin polarized currents. For instance, 
a double dot structure driven by ac-ficld in the presence 
of magnetic field turn out to be a robust spin filtering and 
pumping device.—"^ By applying oscillating gates (radio 
frequency) in an open quantum dot in the presence of 
Zeeman field, an adiabatic spin pump was gcneratedi^"— 



The current ringing^ that arises in a quantum dot system 
when a bias voltage is suddenly switched on develops spin 
dependent beats when the dot level is Zeeman split i^i^ 
Coherent quantum beats in the current, spin current and 
tunnel magnetoresistance of two dots coupled to three 
ferromagnetic leads were also reported recentlyi^ Addi- 
tionally, spin spikes take place when a bias voltage is 
abruptly turned off in a system of a QD attached to fer- 
romagnetic leads 

The study of quantum transport of spin polarized elec- 
trons in the presence of time varying fields was greatly 
motivated by the development of experimental tech- 
niques. These techniques allow for the coherent control of 
the complete dynamics (initialization-manipulation-read 
out) of single electron spins in quantum dotsi^ Particu- 
larly, some of those used to coherently manipulate spin 
states are based on time dependent gate voltages<^"— 

In the present work we consider a single level quan- 
tum dot coupled to left and right ferromagnetic leads in 
the presence of a static bias voltage and sinusoidal gate 
voltage. The oscillating gate potential introduces addi- 
tional photon-assisted conduction channels that can be 
tunned via a dc-gatc field to lie within the conduction 
window of the system. In the presence of Zeeman split- 
ting, produced by an applied magnetic field, the contribu- 
tion from the photon-assisted channels becomes different 
for spins up and down, resulting in photon-assisted spin 
polarized currents. It is worth mentioning that this effect 
takes place even in the absence of ferromagnetic leads. 
However, when the leads are ferromagnetic and paral- 
lel aligned, the resonant current peaks are amplified for 
one spin component and suppressed for the other. Thus 
the photon-assisted current-polarization is enhanced. We 
also calculate the tunnel magnetoresistance (TMR) as a 
function of the gate frequency, which exhibits a variety of 
peaks and dips, having even a changed of sign, depending 
on the gate frequency. 

The paper is organized as follows: in Sec.|ll]we present 
the theoretical model and describe the formulation based 
on nonequilibrium Green's function technique and in 
Sec. mil we show and discuss the numerical results. Fi- 
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FIG. 1. (color online) Energy diagram for the system con- 
sidered. A quantum dot is coupled to a left and to a right 
ferromagnetic electron reservoirs via tunneling barriers. The 
dot level is Zeeman split. A capacitively coupled gate terminal 
introduces a periodic perturbation of the dot level. This mod- 
ulation induces additional photon-assisted channels (dashed 
lines) for spin polarized transport. 



nally, in Sec. IIVI we present our concluding remarks. 



II. MODEL AND THEORETICAL 
FORMULATION 

For concreteness, the energy profile of our system is 
illustrated in Fig. [T] and is described by the Hamiltonian, 
H = Hl+ Hr + Hoit) + Ht, where 

describes the free electrons in the (L) or the right (R) 
lead, in which CkcrL(i?) [Cko.L(_R)] ^® operator that an- 
nihilates [creates] an electrons in the lead L{R) with mo- 
mentum k, spin a and energy ekCTL(_R)- We consider a 
static source-drain applied voltage (eVsn ~ I^l — M-r) 
which drives the system out of equilibrium, breaking the 
left/right symmetry of the Hamiltonian. The time de- 
pendence of our Hamiltonian is fully accounted via the 
dot Hamiltonian, 

Hz5W = 5^e.W4da, (2) 

a 

where e^{t) = ed{t) + aEz/2, with ed{t) being the time- 
dependent dot level and Ez a Zeeman splitting of the dot 
level due to an external magnetic field. Here we use a = 
-|- and (T = — for spins up and down, respectively. The 



operator (dj-) annihilates (creates) one electron with 
spin a and energy ea-{t) in the dot. In practice the time 
dependence in the dot level is controlled by an oscillating 
gate voltage Vg{t), such that ed{t) = eo + eVg{t), where cq 
is the dc component of the energy and eVg{t) = Acos(wi) 
oscillates with amplitude A and frequency uj. Finally 

ijT = ^(v^4/. + i^*4ck.„), (3) 

describes the tunnel coupling between the leads and the 
dot, with a constant coupling strength V and allows for 
current to flow across the QD. 

To calculate the time dependent spin polarized current 
we employ the Keldysh Green's function formalism^^ that 
allows for an appropriate approach to our noncquilibrium 
time-dependent situation. Starting from the current def- 
inition 12{t) = -e{N„) = -ie{[H,N„]}, where N„ is the 
total number of particle operator for spin cr (here we take 
h= 1), the current can be written as^ 

J^(<)=2eRe|5]FG<,,^(i,0|, (4) 

where G< j,^^^(i,t') = i{c\^^{t')d„{t)) . Using the equation 
of motion technique and taking analytical continuation^ 

to obtain G^, (t, t') one finds to the current the follow- 
er, ^ctt/ ^ ' ' 

ing 

T2{t) = -2er;im{ j ^ f Jhe-^<'^-'^ 

X [Gl,{tM)U^) + G<,{tM)]}. (5) 

where /r;(e) is the Fermi distribution function of the 7;-th 
lead, and TJJ = 2'k\V\'^ gives the tunneling rate between 
lead rj and dot for spin component a. is the density 
of states for spin a in lead f]. In the present model we 
assume constant density of states (wide-band limit). The 
ferromagnetism of the electrodes is modeled by consid- 
ering rjj = ro(l ±p^) where -I- (— ) stands for spin up 
(down), is the polarization of lead yy-t h^^i'^^ and Fq the 
tunneling rate strength. The quantity Fq is fixed along 
the paper, so all the other energies will be expressed in 
units of Fq. We consider both parallel (P) and antipar- 
allel (AP) alignments of the lead polarizations. In the P 
case we assume majority down population in both leads, 
while in the AP configuration we take majority down 
population in the left lead and majority up population 
in the right lead. In terms of the parameters we have 
PL ~ PR = P = —0.4 for the P and pl = —pR = p = —0.4 
for the AP case4S Taking the time average of the current 
we find 

(6) 

where 

°° A 

{A,{e,t))= Jn{-)9t(^,u;), (7) 
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FIG. 2. (color online) Color map of the total transmission 
coefficient T{e,ij) = T-|-(£, tj) + Tj,(e, a;) as a function of fre- 
quency and energy in the parallel alignment (pL = Pr = 
—0.4). For increasing lj, T{e,uj) develops additional photon- 
assisted peaks that allows off-resonant spin transport. The 
main two central peaks correspond to the Zeeman split levels 
= eo + Ez/2 and — eo — Ezl2. The satellite peaks are 

given by e|"' = eo + Ez/2 ± nuj and e|"' = eo — Ez/2 ± nu), 
with n = 1, 2, 3, .... For increasing uj, the satellite peaks tend 
to vanish and the system recovers its original two levels and 
ej. The horizontal dashed lines delimit the conduction win- 
dow [p.r,^l\- Units: Energy in units of To and ujq = To/h. 
Parameters: eo = — ISFo, Ez — 4ro, A — SFo, PL = IFo, 
^J,R = 0. 



with J„ being the n-th order Bessel function and 
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] ^ . Here we used the 



-I- Acos(wt), with 6° = eo + a Ez/2. 
Substituting this result into Eq. the current can 
be written in its Landauer form^ 

dt. 
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FIG. 3. (color online) Multiplet structure developed in the 
presence of an oscillating gate frequency. The black lines cor- 
respond to spin up while the gray lines to spin down. The 
levels are shifted linearly with the gate frequency, following 
eS"' ~ e° zbntj, n = 1, 2, 3, ... The up and down levels are Zee- 
man split. The horizontal dashed lines correspond to the left 
(/ii) and to the right (/ia) chemical potentials. The channels 
e|"' and ej"' attain resonance within the conduction window 
[m^i Mfl] foi" certain frequencies, which differ for each spin com- 
ponent. Units: Energy levels in units of Fq and loq = Fo/fi. 
Parameters: eo = — 15Fo, Ez = 4Fo, A = 5Fo, PL = IFo, 
Atfl = 0. 



functions are approximated by step functions. In this 
regime, the integral in Eq. is carried out in the range 
[^i , Mi?] , thus resulting in 



(10) 



where = e pL°^pH is the resonant current without mod- 
ulated gate voltage and 

E ^,'.(--)[e,f H-e,^«(c.)]/7r, (11) 



n— — oo 



with Q"'{uj) ~ arctan[2(/.t,, — e° — nLij)/Tcr]- In what 
follows we present our numerical results to the spin po- 
larized transport. 



III. NUMERICAL RESULTS 



where e^" = e" + mu. Eq. ^ shows that the harmonic 
modulation of the dot level yields to photon-assisted 
peaks in the transmission coefScient.— In addition to this, 
here we have the spin splitting of these peaks and the 
ferromagnetic leads, that results in an enhanced spin 
photon-assisted transport. 

A further simplification can be made in Eq. ([5]) by 
considering the low temperature regime, where the Fermi 



Figure [2] shows the sum T = T|- Tj_ as a function 
of uj and energy in the case of polarized leads with par- 
allel magnetizations. As ui increases, a multiplet struc- 
ture takes place in the transmission coefficient [Eq. 
The two central peaks in r(e, uj) correspond to and ej, 
while the lateral peaks are related to e° ± nuj. Due to 
the Zeeman splitting, the whole pattern for T-]- is shifted 
upward while is moved downward. The highest of the 



peaks are strongly affected by the frequency. For the n-th 
peak its amphtude is given by 4r^rf J^(A/a;)/r^. For 
sufficiently large w, the additional photon-assisted peaks 
are suppressed, remaining only the two central peaks. 
The broadening difference for up and down spin channels 
comes from the ferromagnetism of the electrodes that are 
parallel aligned, with majority down population in both 
sides [pl ~ Pr — —0.4). This gives rise to narrower spin 
up peaks than spin down ones. In the case of parallel 
aligned leads with majority up populations, we have ba- 
sically the same structure, but with an inversion of the 
peak widths, with spin up peaks now becoming broader. 

Spin polarized transport can arise depending upon the 
position of the peaks of ^^(e, oj) with respect to the con- 
duction window. Fig. ([3]) shows the channels e° -I- nuj 
and e° -I- nui for n = 0,±1,±2,±3 as black and gray 
lines, respectively. The left and right chemical potentials 
are indicated by the horizontal dashed lines. A net elec- 
tron transport from the left to the right lead can take 
place whenever a channel + nuj attains the conduc- 
tion window (CW) interval [fiLifJ-B]- Due to the Zeeman 
splitting, each spin component crosses fi^ or fiR at dif- 
ferent frequencies, thus resulting in a frequency selective 
spin transfer between the leads. In Fig. we indicate 
by up and down arrows the corresponding crossing of the 
CW for spins t and J., respectively. In the present study 
we focus on the off-resonant regime, where the dot levels 
e!j! and ej are below the CW. In this case only photon- 
assisted electrons can tunnel through the system. In or- 
der to match this condition we adopt to the numerical 
parameters the following values: eq ~ — ISFq, Ez = 4ro, 
A = SFq, /ii = IFq and = 0. Later on we will 
also look at distinct parameters in order to explore the 
robustness of our main results. In experiments we find 
typically Fq ~ lOOneVM^— So to the parameters as- 
sumed we have Ez ~ AOOneV. This Zeeman energy 
split is reasonable for semiconductor quantum dots in the 
presence of magnetic fields ^ 1-10 T4^ Additionally, for 
these values we find wg = ^ ~ 150 GHz. So the present 
theoretical effects could be observed for gate frequencies 
around 1.5 THz [oj ~ lOwo, see Fig. (|4])].ii Alterna- 
tively, if Fq is reduced to ^ l/^eV,— we obtain gate 
frequencies around w ~ IOwq ~ 15 GHz, which is quite 
feasible experimentally.— Our currents will be given in 
units of Iq = eVo/h, which is in the range Iq ^ 0.24nA- 
24nA for Fq ^ IficY - 100 ficV. Since our spin resolved 
photon-assisted currents are typically 10"'^ /q, wc have 
pA currents, which could be measured with picoampere 
measurement technologies. 

Comparing Fig. ([3|) to Fig. ([2]) one can note that 
even though n = 3 and n = 2 attain resonance with- 
ing [fj,L,^fi], their corresponding transmission amplitude 
are very low, which makes the transport weak via those 
channels. In contrast, the n = I channels, for up and 
down spins, have a higher transmission amplitude, which 
makes the spin transfer via these channels more appre- 
ciable. 

Fig. (m shows the up and down components of the 
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FIG. 4. (color online) Spin resolved currents against gate fre- 
quency for leads (a) nonmagnetic and (b)-(c) ferromagnetic 
(black lines for spin up and gray lines for spin down). In 
panels (b) and (c) we show the parallel and antiparallel align- 
ments, respectively. Both up and down currents show peaks 
corresponding to the crossing of e|"^ and e|"^ illustrated in 
Fig. ((3]). The highest peak for each spin component comes 



from the resonance of the levels ej^' and e^'' , withing the con- 
duction window. In the parallel alignment the majority down 
population in both leads turns into an amplification of the 
down current. In the antiparallel case, though, the currents 
are very similar to the nonmagnetic case. Units: Iq = eFo/fi 
and iJo — To/h. Parameters: eo = — ISFo, Ez = 4ro, 
A = SFo, fiL = IFo, Mfl = 0, p = -0.4. 



current against gate frequency. Three cases are consid- 
ered: (a) nonmagnetic leads, ferromagnetic leads in the 
(b) parallel and (c) antiparallel alignments. In all the 
three cases two major peaks are found {n = 1). Satel- 
lite peaks for high order channels {n ~ 2, 3) are also 
seen. Each peak emerges whenever a channel + nuj 
enters the conduction window [indicated by f and J, ar- 
rows in Fig. ([3])]. Due to the Zeeman splitting, the res- 



onance for spin up arises in lower frequencies than that 
for spin down. Additionally, the interplay between Zee- 
man splitting and the amplitude of the transmission co- 
efficient results into a higher peak for spin up than for 
spin down in the case of nonmagnetic leads?^ Further 
amplification of the spin down current peak is observed 
when the leads are made ferromagnetic. In Fig. Sl^b) 
we present and for leads parallel aligned, with a 
majority down population in both sides. This means 
that we assume for the polarization parameter a nega- 
tive value, with = = p = —0.4. This implies that 
pL,fl > F:^'^ iPi'^ > which favors more the spin 

down electrons to tunnel through the system, thus in- 
creasing the spin down current peak. In the antiparallel 
case, where we have a majority down population in the 
left lead and a majority up population in the right lead 
(T^ > F^ and Tf < for pL = -pR = P = -0.4), 
the incoming and outgoing rates compensate each other. 
This results into equal weights for both up and down cur- 
rents, so the current remains almost the same compared 
to the nonmagnetic case. 

Fig. ([5]) shows how the spin resolved currents evolve 
when the leads polarizations are enlarged in both P and 
AP configurations. When p becomes more negative, the 
spin down tunneling rates and are strengthened, 
while and are diminished in the parallel case. 
This amplifies the peak for spin down current while sup- 
presses the peak for spin up, thus making the current 
more down polarized. Eventually, for larger enough \p\ 
the current dominates over /-[■ for all gate frequen- 
cies. Conversely, in the antiparallel configuration, as p be- 
comes more negative, both and /j. are suppressed. One 
may note that when we pass from P to AP alignment, 
the major spin up peak has its width broadened. This 
can also be seen by calculating the width of this peak: 
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for P and AP, respectively. In particular, this broadening 
effect makes the total antiparallel current slightly higher 
than the parallel current (/^ < I^^), and this results 
into a negative magnetoresistance, as we will see next. 

In Fig. (|6]) we show the total current in both P and 
AP alignment and the tunnel magnetoresistance, defined 
according to 



TMR = 



I^-I 



AP 



jAP 



(12) 



where I^/^p = I^^^^ + I^^^^ . The highest current 

peak is related to the resonance of the channel e^"^ with 
the CW. The second highest peak is due to spin down 



resonance, e^^'^ ~ fip. The spin valve effect is clearly 
seen for almost all frequencies, i.e., > I^^ . However, 
for frequencies around lOwo, we observe < I^^ for 
p = —0.4 [Fig. Elja)]. This photon-assisted opposite spin 
valve effect is related to the broadening of the spin up 
peak discussed in Fig. ([5]), when the lead polarizations 
rotate from P to AP alignment. This feature is reflected 
in the TMR, which acquires negative values. Notice that 
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FIG. 5. (color online) Spin resolved currents against gate 
frequency for different polarization p, in both parallel and 
antiparallel configurations (black lines for spin up and gray 
lines for spin down). In the plots p is negative, which means 
that both leads have majority spin down population in the 
P case and majority down (up) population in the left (right) 
lead for the AP case. In the P alignment when p increases (in 
modulus) the tunneling rates between the dot and the leads 
enlarge for spin down and reduce for spin up. This results 
in an amplification of /j. and a suppression of as observed. 
In the AP alignment both and are suppressed as \p\ 
increases. Units: 7o = eTo/h and cjq — To/h. Parameters: 
eo = — ISFo, Ez = 4ro, A = SFq, fJ-L = IFq, = 0. 



near the position of the spin up peak the TMR is fully 
suppressed, while near to the peak of the spin down cur- 
rent it is enlarged for all values of p. This indicates that 
the magnetoresistance is mainly dominated by spin down 
transport. 

Since from the experimental point of view the quan- 
tities VsD, and Ez are relative easily tunable ex- 
perimental parameters, we analyze how these quantities 
affect the present results, aiming to providing a better 
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FIG. 6. Total current + I^) in both P and AP alignment 
and TMR as a function of the gate frequency. From panels 
(a) to (c) the left and right leads polarization are enlarged 
(p = —0.4, —0.6, —0.8). The highest current peak is pre- 
dominantly due to spin up transport while the second highest 
peak is more spin down like. The TMR tends to be sup- 
pressed or amplified around the spin up or spin down peaks, 
respectively. In particular, the TMR attains negative val- 
ues for p = —0.4 in a short gate frequency range, due to a 
photon-assisted inversion of the spin valve effect (/^ < r^^). 
Units: Jo = eVo/h and ujq — Fo/fi. Parameters: eo = — 15Fo, 
Ez = 4Fo, A = 5Fo, /il = IFo, m« = 0. 



FIG. 7. (color online) Spin resolved currents against gate fre- 
quency for both parallel and antiparallel configurations (black 
lines for spin up and gray lines for spin down). From panels 
(a) to (c) we increase the bias voltage [(a) eVsD = SFq, (b) 
eVsD = lOFo, (c) eVsD = ISFo], thus enlarging the conduc- 
tion window [dashed lines in Fig. [3]. Fig. EJa) is similar to 
Fig. [5ja), except by the width of the peaks that are slightly 
enlarged by the CW. When eVs d increases even further [pan- 
els (b)-(c)] the peaks turn even more broaden. In the insets 
we show the TMR for each bias voltage. For eVsD = 5Fo 
the TMR presents a negative value around uj « 12.5tJo- For 
eVsD = lOFo and eVsD = ISFo the TMR becomes positive 
for all frequencies. Units: /o = eFo/fi and luq = Fo/h. Pa- 
rameters: eo = — ISFo, Ez = 4Fo, A = SFq, fiL—tJ,R = eVso, 
p = -0.4. 



guide for future experimental realizations. For instance, 
using the same set of parameters adopted previously 
[pL = -0.4, PR = -0.4(P), pH = +0.4(AP), eo = -WFo, 
A ~ SFo, Ez = 4ro], in Fig. ^ we show how the spin 
resolved currents evolve when the source-drain voltage 
increases in both parallel and antiparallel configurations. 
The bias voltages considered are (a) eVsD = SFq, (b) 
lOFo and (c) ISFq. These values were also indicated by 
dashed lines in Fig. By comparing Fig. ^a) to Fig. 



[5ja) one can note that the peaks are broadened as Vsd 
increases. This feature can be understood looking at the 
different conduction windows (dashed lines) in Fig. 
As fj.L — ^^R. becomes larger, the frequency range in which 
the photon-assisted channels remains inside the CW be- 
comes wider. In the insets of Fig. ([7]) we show the TMR 
against frequency for each Vsd considered. Notice that 
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for eVsD = IOFq and eVso = ISFq [panels (b) and (c), 
respectively] the TMR is positive for all frequencies. 

In order to obtain a clear picture on the set of pa- 
rameters needed to obtain a negative TMR, in Fig. (|8]) 
we plot in a color map of the total transmission coeffi- 
cient difference, AT = — T^^, against gate frequency 
(horizontal axis) and energy e (vertical axis) , for two dif- 
ferent Zeeman splitting energies: (a) Ez = 4Fo and (b) 
Ez = STq. Notice the appearance of negative values 
(dark regions) around the spin up peaks, when the up 
and down branches are not overlapping. The overlap can 
be avoided by controlling the Zeeman splitting. Observe 
that for Ez = SFq [panel (b)] the spin up and spin down 
branches are far enough from each other to ensure rela- 
tively large dark areas. This effect can be understood in 
terms of the broadening of the photon-assisted peaks. In 
the P configuration, the spin down peaks at the transmis- 
sion coefficient are broader than the spin up ones, as can 
be seen by the rates F^ = F;['-hFf = 2Fo(l + |p|) = 2.8Fo 
and F^ = 2Fo(l - |p|) = 1.2Fo {\p\ = 0.4). So when 
these peaks coincide at the same energy {Ez = 0) the 
spin down peaks lie above the spin up ones, thus mak- 
ing the total transmission coefficient dominated by the 
spin down component. When the magnetic alignment is 
rotated from P to AP, the width of the spin up peak in- 
creases and spin down diminishes, becoming both 2Fo. 
This facilitates T^^ > and Tf^ < nearby each 
transmission peak. If the peaks are far enough from each 
other (Ez ^ 0) it is possible to obtain T = Tt + « Tt 
around the up peaks, leading to T^^ « T^^ > « 
[see Fig. M,c) for clarity] . This is the main condition for 
a negative TMRi^ For instance, by tunning the set of 
parameters such that the CW lies on the darker area of 
the map, the TMR becomes negative. By increasing the 
CW, we can eventually cover an energy range with more 
bright than dark areas, resulting in a positive TMR. 

In Fig. (jni we show the current vs. uj obtained for 
the conduction window set as in Fig. |Hl[b) [dashed (green) 
lines] Since the CW is dominated by a dark region 
(T^ < T^P) we expect TMR < 0. The spin up and spin 
down currents in Fig. |9l[a) reveal contrasting behavior. 
While the spin down current shows the same behavior 
already seen in previous results [e.g.. Fig. ([S])], the spin 
up current oscillates and then it increases to a higher 
saturation value. This higher value of I-^ in the large fre- 
quency limit can be understood by looking at the trans- 
mission coefficient in panels EJb) and (c) . Spin up and 
spin down transmissions coefficients are drawn in both P 
and AP configurations. The vertical thin solid lines give 
the border of the CW. Comparing the amplitude of 
inside the CW for both lj <^ ljq and a; 3> wq, we observe 
that T-^ becomes amplified for larger frequencies, which 
results into higher spin up current in this limit. Con- 
versely, the spin down transmission coefhcient is slightly 
higher in the CW for a; <C wq, resulting in a spin down 
current a bit higher in the low frequency limit, compared 
to its value in the high frequency regime. 

To better understand the oscillatory structure found 
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FIG. 8. (color online) Two-dimensional map of the difference 
AT = — against frequency a; and e for two different 
Zeeman splitting: (a) Ez — 4ro and (b) Ez — SFq. The dark 
regions indicate > which results in negative TMR. 

As Ez increases [from panel (a) to (b)] the dark ares enlarge. 
This is so because the spin up and down channels become 
apart from each other, allowing to have = T:^'' + T^^ ~ 
T^^ > = + rf « Tf (broadening effect) around the 
spin up channels. By playing with the set of parameters one 
can place the conduction window of the system in the darker 
region of the map, which results into a more negative TMR. 
Units: Iq = eVo/h and ij;o = Fo/fi. Parameters: eo = — 4.5ro, 
A = SFu, /ii = IFo, Mfl = 0, p = -0.4. 



in the spin up current, we plot in Fig. the spin up 
parallel transmission coefhcient vs. energy for various 
cj. As cj increases the transmission coefficient develops 
a variety of peaks mainly in the range between e!^ — 5Fo 
and 6° -I- 5Fo (e° ± A) . Since this energy interval contains 
the CW, some of the photon-assisted peaks that arise for 
increasing w appear within the CW and eventually lie 
outside it for large enough followed by a new peak 
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FIG. 9. (color online) (a) Spin resolved currents in both P 
and AP alignments against gate frequency for the conduction 
window drawn in Fig. [8jb). Similar to previous figures, the 
spin down component reveals peaks at some particular fre- 
quencies due to the matches ei"^ = fin. In contrast, the spin 
up current oscillates and then it tends to a relatively large sat- 
uration value for large uj. In the inset of Fig. [9ja) we show the 
TMR which acquires negative values in accordance to the CW 
in Fig. [HJb). Panels (b)-(c) present the spin resolved trans- 
mission coefficient for (b) small and (c) large u in both P and 
AP configurations. The thin vertical lines denote the border 
of the CW. In particular, observe in panel (c) that inside the 
CW we have T^^ > . Since T^^"^^ is too small at the CW, 
the total transmission coefficient T^^^^ is given essentially by 



rf/^^,soT- 



around the peak. This results in a negative TMR. In pan- 



els (d)-(e) we plot (d) Tf 



and (e) for different values 



of Lo. The curves were vertically displaced for clarity. As 
u) increases additional peaks emerge mainly in the interval 
± A = e° ± SPo . More specifically, in the spin up case the 
appearance and suppression of these additional peaks inside 
the CW gives rise to the oscillatory patter of in the low 
frequency limit. Units: Jo = eVo/h and ujq = Po/ft. Parame- 
ters: eo = -4.5r,j, Ez = SFo, A = SPo, fiL = IPo, ^J^.R = 0, 
p = -0.4. 



emerging inside it. This results in an oscillatory pattern 
of the current. In panel |9l[c) we show as function of 
energy for the same values of uj as in panel (d). Since 
ej + A < fifi the photon-assisted peaks can only cross 
the CW for particular frequencies, which gives rise to 
the peaks observed. Finally, in the inset of Fig. [9l[a) we 
plot the TMR, which presents very low negative values 
-30%). This negative TMR is a result of the set of 
parameters chosen from Fig. |8l 



IV. CONCLUSION 



We have studied spin polarized transport in a quantum 
dot attached to ferromagnetic leads in the presence of an 
osciUating gate voltage Vg{t). A static source-drain bias 
voltage is also applied in order to generate current. The 
oscillating Vg{t) gives rise to photon-assisted transport 
channels that allow electrons to flow through the system. 
Due to a Zeeman splitting of the dot level, the photon- 
assisted contributions to the transport are distinct for 
spins up and down, providing an interesting way to ob- 
tain current polarization that can be controlled by gate 
frequency. As the leads polarization is enlarged, with a 
majority down population in both leads (P alignment), 
the spin down photon-assisted current peak is enhanced, 
while the spin up peak is suppressed. Moreover, when the 
relative polarization alignment of the leads is switched 
from P to AP, the width of the main spin up peak of the 
current is broadened. This additional broadening effect 
results in an opposite spin valve behavior (^I^ < I^^) 
for gate frequencies around the spin up resonance. As 
a result, a photon-assisted negative tunnel magnetoresis- 
tance is found. 



ACKNOWLEDGMENTS 



The authors acknowledge A. P. Jauho and J. M. Villas- 
Boas for valuable comments and suggestions. This work 
was supported by the Brazilian agencies CNPq, CAPES 
and FAPEMIG. 



^ N. S. Wingreen, A. P. Jauho, and Y. Meir, Phys. Rev. B ^ L. J. Geerligs, V. F. Anderegg, P. A. M. Holweg, J. E. 
48, 8487 (1993). Mooij, H. Pothier, D. Esteve, C. Urbina, M. H. Devoret, 



9 



Phys. Rev. Lett. 64, 2691 (1990). 
^ L. P. Kouwenlioven, A. T. Johnson, N. C. van der Vaart, 
C. J. P. M. Harmans, C. T. Foxon, Phys. Rev. Lett. 67, 
1626 (1991). 

* C. Bruder and H. Schoeller, Phys. Rev. Lett. 72, 1076 
(1994). 

^ L. P. Kouwenhoven, S. Jauhar, K. McCormick, D. Dixon, 
P. L. McEuen, Yu. V. Nazarov, N. C. van der Vaart, and 

C. T. Foxon, Phys. Rev. B 50, 2019(R) (1994). 

® L. P. Kouwenhoven, S. Jauhar, J. Orenstein, P. L. McEuen, 
Y. Nagamune, J. Motohisa, and H. Sakaki, Phys. Rev. 
Lett. 73, 3443 (1994). 

B. J. Keay, S. J. Allen, Jr., J. Galan, J. P. Kaminski, K. 
L. Campman, A. C. Gossard, U. Bhattacharya, and M. J. 
W. Rodwell, Phys. Rev. Lett. 75, 4098 (1995). 

* For a review see G. Platero and R. Aguado, Phys. Rep. 
395, 1 (2004). 

® C. A. Staflord and N. S. Wingreen, Phys. Rev. Lett. 76, 
1916 (1996). 

1° A. F. Amin, G. Q. Li, A. H. Philhps, and U. 
Kleinekathofer, Eur. Phys. J. B 68, 103 (2009). 
F. Cavaliere, M. Governale, and J. Konig, Phys. Rev. Lett. 
103, 136801 (2009). 

R. Zhu and H. Chen, Appl. Phys. Lett. 95, 122111 (2009). 
E. Prada, P. San- Jose, and H. Schomerus, Phys. Rev. B 
80, 245414 (2009). 
" G. M. M. Wakker and M. Blaauboer, Phys. Rev. B 82, 
205432 (2010). 

M. Alos-Palop and M. Blaauboer, arXiv: 1102.0926 (2011). 
J. Splettstoesser, M. Governale, J. Konig, and M. Biittiker, 
Phys. Rev. B 81, 165318 (2010). 

S. Kurth, G. Stefanucci, E. Khosravi, C. Verdozzi, and E. 
K. U. Gross, Phys. Rev. Lett. 104, 236801 (2010). 
^® G. A. Prinz, Science 282, 1660 (1998); S. A. Wolf, D. 

D. Awschalom, R. A. Buhrman, J. M. Daughton, S. von 
Molnar, M. L. Roukes, A. Y. Chtchelkanova, and D. M. 
Treger, ihid. 294, 1488 (2001); Semiconductor Spintron- 
ics and Quantum Computation, edited by D. Awschalom, 

D. Loss, and N. Samarth (Springer, Berlin 2002); D. D. 
Awschalom and M. E. Flatte, Nat. Phys. 3, 153 (2007). 

" E. Cota, R. Aguado, and G. Platero, Phys. Rev. Lett. 94, 
107202 (2005). 

R. Sanchez, E. Cota, R. Aguado, and G. Platero, Phys. 
Rev. B 74, 035326 (2006); R. Sanchez, E. Cota, R. Aguado, 
and G. Platero, Physica E 34, 405 (2006). 

E. R. Mucciolo, C. Chamon, and C. M. Marcus, Phys. Rev. 
Lett. 89, 146802 (2002). 

2^ T. Aono, Phys. Rev. B 67 155303 (2003). 

S. K. Watson, R. M. Potok, C. M. Marcus, and V. Uman- 
sky, Phys. Rev. Lett. 91, 258301 (2003). 

F. M. Souza, Phys. Rev. B 76, 205315 (2007). 

E. Perfetto, G. Stefanucci, and M. Cini, Phys. Rev. B 78, 
155301 (2008). 

P. Trocha, Phys. Rev. B 82, 115320 (2010). 

F. M. Souza, S. A. Leao, R. M. Gester, and A. P. Jauho, 
Phys. Rev. B 76, 125318 (2007). 

F. M. Souza and J. A. Gomez, Phys. Status Sohd B 246, 
431 (2009). 

R. Hanson and D. D. Awschalom, Nature 453, 1043 
(2008). 

J. M. Elzerman, R. Hanson, L. H. W. van Beveren, B. 
Witkamp, L. M. K. Vandersypen, and L. P. Kouwenhoven, 
Nature 430, 431 (2004). 



R. Hanson, L. H. W. van Beveren, I. T. Vink, J. M. Elz- 
erman, W. J. M. Naber, F. H. L. Koppens, L. P. Kouwen- 
hoven, and L. M. K. Vandersypen, Phys. Rev. Lett. 94, 
196802 (2005). 

A. C. Johnson, J. R. Petta, J. M. Taylor, A. Yacoby, M. D. 
Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Nature 435 925 (2005). 

J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 
Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, A. C. 
Gossard, Science 309, 2180 (2005). 

F. H. L. Koppens, C. Buizert, K. J. Tielrooij, L T. Vink, 
K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and L. M. 
K. Vandersypen, Nature 442, 766 (2006). 
C. Barthel, J. Medford, C. M. Marcus, M. P. Hanson, and 
A. C. Gossard, Phys. Rev. Lett. 105, 266808 (2010). 
^® H. Haug and A. P. Jauho, Quantum Kinetics in Trans- 
port and Optics of Semiconductors, Springer Solid-State 
Sciences 123 (1996). 

In the present calculation we follow the formulation devel- 
oped in A. P. Jauho, N. S. Wingreen, and Y. Meir, Phys. 
Rev. B 50, 5528 (1994). 

W. Rudzihski and J. Barnas, Phys. Rev. B 64, 85318 
(2001). 

F. M. Souza, J. C. Egues, and A. P. Jauho, Braz. J. Phys. 
34, 565 (2004). 

In the present work we consider only the P and AP align- 
ments. We expect that for more general non-coUinear mag- 
netizations the present currents will vary monotonically 
from their parallel (0 = 0) to their antiparallel (0 = tt) val- 
ues. In the adiabatic-pumping regime it was found (see, for 
instance, Ref. [4l[) a monotonic suppression of the pumped 
charge as the angle between the magnetizations of the left 
and right leads is continuously varied from = to </) = tt. 
J. Splettstoesser, M. Governale, and J. Konig, Phys. Rev. 
B 77 195320 (2008). 

^"^ An alternative but equivalent version of this equation was 
derived in Ref. , where the distinction between this result 
and the Tien-Gordon formula was discussed. 

^'^ D. G.-Gordon, H. Shtrikman, D. Mahalu, D. A.-Magder, 
U. Meirav, M. A. Kastner, Nature 391, 156 (1998). 

^'^ D. Goldhaber-Gordon, J. Gores, M. A. Kastner, H. Shtrik- 
man, D. Mahalu, and U. Meirav, Phys. Rev. Lett. 81, 5225 
(1998). 

F. Simmel, R. H. Blick, J. P. Kotthaus, W. Wegscheider, 
and M. Bichler, Phys. Rev. Lett. 83, 804 (1999). 

*® R. Hanson, B. Witkamp, L. M. K. Vandersypen, L. H. 
Willems van Beveren, J. M. Elzerman, and L. P. Kouwen- 
hoven, Phys. Rev. Lett. 91, 196802 (2003). 

"^"^ Significant improvements on THz electronics have been re- 
cently achieved. See for instance, T. W. Crowe, W. L. 
Bishop, D. W. Porterfield, J. L. Hesler, R. M. Weikle, IEEE 
J. Sohd-State Circuits, 40 2104 (2005). 
S. Amasha (private communication). 

''''' K. MacLean, S. Amasha, I. P. Radu, D. M. Zumbiihl, M. 
A. Kastner, M. P. Hanson, and A. C. Gossard, Phys. Rev. 
Lett. 98, 036802 (2007). 

5° H. Eisele, S. P. Khanna, and E. H. Linfield, Appl. Phys. 
Lett. 96 072101 (2010). 

Two split peaks in a pumping current, one peak for each 
spin component, were also found by Cota et al. in Ref. [l^ 
for a double dot system in the presence of magnetic field. 
In the case of majority up population in both leads in the 
P alignment, we have the opposite behavior seen in Fig. 
((8]), i.e., the maxima become minima and vice- versa. So 



10 

the dark areas in the plot (T^ < T^^) arise around the spin down branches, instead of the spin up ones. 



